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Q\ • Abstract 

The zero modes and current carrying capability of a cosmic string 
formed at one phase transition can be modified at subsequent phase 
O || transitions. A new, generalised index theorem is derived that is ap- 

plicable to theories with multiple phase transitions. This enables us 
to investigate the fate of string zero modes during sequences of phase 
transitions in a variety of models. Depending on the couplings that 
the breaking introduces, the zero modes may be destroyed and the 
superconductivity of the string removed, and thus vortons dissipate. 
We discuss the features of the theory that are required to produce this 
behaviour and consider the implications of spectral flow. 
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1 Introduction 



Topological defects, and in particular cosmic strings, arise in many grand 
unified theories They are likely to be copiously produced during phase 
transitions in the early universe. A network of cosmic strings could explain 
the observed anisotropy in the microwave background radiation and the 
large scale structure of the Universe In the past few years it has been 
realised that they may have considerably richer microstructure than previ- 
ously realised M. In particular, the presence of fermion zero modes || in 
the spectrum of a cosmic string has profound implications for the cosmology 
of the defects. Zero modes render the string superconducting [f|, allowing 
currents of up to 10 20 Amps to flow along the string. Stable relics, vortons, 
can form as collapsing string loops are stabilised by the angular momentum 
of the trapped charge carriers ||. These relics can overdose the Universe 
unless stringent constraints are applied to the string model Q . 

The cosmological implications of the vortons are most pronounced when 
the Universe has become matter dominated. If they decay during the era of 
radiation domination the cosmological catastrophe may be avoided and the 
vorton bounds evaded. However, it has recently been realised that subse- 
quent phase transitions can have a considerable effect on the microphysics 
of cosmic strings. Zero modes on cosmic strings can be both created and 
destroyed ||, thus creating or destroying currents on strings. In this paper 
we discuss the fate of vortons and string superconductivity as the strings 
encounter subsequent phase transitions in a systematic fashion. 

In section |2] we derive an index theorem giving the number of zero modes 
for a general mass matrix. Whilst index theorems have been derived before 
|]. [l(J they have been more restrictive in their validity and have only been 
able to determine the difference in right-moving and left-moving zero modes. 
Our index theorem has much more general applicability and can give a bound 
on the number of zero modes. 

In section [3] we discuss the fate of zero modes on strings formed at the 
breaking of an SO(10) Grand Unified symmetry ||. At the electroweak 
phase transition these zero modes acquire a small mass which leads to dissi- 
pation of the string current. This allows vortons to decay and weakens the 
cosmological bounds on such models [11]. 

In other models the zero modes survive subsequent transitions allowing 
the associated vortons to persist. Such behaviour is displayed by a toy model 
discussed in section f|, where we explicitly construct the zero mode solutions 
of the Dirac equation after the symmetry breaking. 
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We compare these two models in section || and consider the implications 
of spectral flow. An important feature that allows zero modes to be removed 
is the presence of a particle that mixes with the zero mode after the tran- 
sition. The implications of such couplings for current build up before the 
transition are also considered. Finally, we summarise our conclusions. 

2 Zero Mode Index 

Cosmic strings form in models with vacuum manifolds which are not simply 
connected. For example in a £7(1) model, with potential (|</>| 2 — t] 2 ) 2 , stable 
solutions exist with <j) = r\e %e at r = oo. In order for the total energy 
to be finite, a non-zero gauge field is needed to give a vanishing covariant 
derivative at r = oo. In a more general theory, involving a larger group, G, 
string solutions take the form 

<Mr,0) = e iT ' e <l>(r) , A e = -T(r) , (1) 

er 

where T s is a generator of G that is broken by <p. The choice of T s is restricted 
by the fact that 4> must be single valued. 0(r) is equal to the usual VEV of <f> 
at r = oo, and must be regular at r = 0. T(r) obeys T(0) = 0, T(oo) = T s . 

In a general theory, T s will affect different components of 4> differently. 
This means that the various parts of <f> can have a wide range of winding 
numbers. In a theory with multiple phase transitions, the additional Higgs 
fields will be affected in the same way. It may also be necessary to alter T s 
at phase transitions to make the new Higgs fields single valued. 

In a theory with nf two-component fermions, the fermionic part of the 
Lagrangian is 

Li = i> a ia»D^ a - -ii, a M aP ^ + (h.c.) (2) 

where t/jg = ia 2 ^J. If M a p depends on 9, as would be expected if M a/ g arose 
from the Higgs field of the string, then it is possible that the field equations 
will have non-trivial zero energy solutions. Solutions with only r and 6 
dependence can be split up into eigenstates of a 3 : ip^, ip^. Such solutions 
have zero energy. If we solve the equations of motion in the background of 
a cosmic string, the field equations become 

e i6 (d r + -d e + eAo) ^ + M a ^f = , (3) 
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e~ ie (d r - l -d e - eA e ^j ^ - M a/ #J* = , (4) 

where Aq is the string gauge field. If z and t dependence is added to the 
solutions they will correspond to currents flowing along the string. Their 
direction is left for those corresponding to (^), and right for (|j). In order to 
be physically relevant the solutions must be normalisable. Let Nl and Nr 
be the number of such solutions to (|3|) and (Q) respectively. We attempt to 
derive an expression for them by generalising the analysis in ref. |3|], which 
involves removing the 9 dependence of the problem, and then considering 
solutions near r = oo and r = 0. 

Choose the ip a s to be eigenstates of the string gauge field, with eigenval- 
ues q a . The q a will depend on the fermion charges and the winding numbers 
of the various components of the Higgs fields. Since the mass terms in (||) 
are gauge invariant, the angular dependence of the mass matrix must be 

M af3 (r,9) = C a p{r)e i{qa+q ? )e (no summation) . (5) 

The 9 dependence Ccin also be factored, out of the ip^s. 

ipa = e l ^-^(U^{r)e ue + Vt{r)e~ iW ) , (6) 
= e ^+\)(jjR{r)e iW + V^*{r)e- lW ) . (7) 

First consider left moving zero modes. Putting (||J7|) into @) gives equa- 
tions for U% and . As r — > oo, C a p = O(l) and eAg = 0(l/r), so 

^ + ^(00)^ = 0, (8) 
d r V£ + C aP (oc)U$ = 0. (9) 

Diagonalising C aj 3 gives 2?if complex solutions proportional to exp(iAo,), 
where X a are C a ps eigenvalues. Thus, assuming all X a 7^ 0, exactly m of 
these large r solutions are normalisable at r = 00. 

The Higgs fields, and hence M a g, are regular at the origin, so as r — > 0, 
C af3 = O(l), and eA e = 0(r). Thus 

+ C af) V£* = 0, (10) 
+ C a/3 U$* = 0. (11) 
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To leading order, the small r solutions are 



where each choice of a = 1 



0(l)r 



V/3, 



(12) 



0{l)r q p-> 2 
nf gives one complex solution, and 



V, 



L 

a 

V L 



0(l)r 
0(l)r 



■i-l 



V/3 , 
V/3 ^ a 



(13) 



which gives a total of 2m independent complex solutions. For given I and a, 



(12) will be normalisable (for small r) if I > — (fa. + 1/2. If Z < (fa — 1/2 then 

(13) will be normalisable. Thus for a given Z the number of well behaved 
small r solutions is 



N° L (l) 



+ 1 



z > 



-9q + 5 



(14) 



where I [X] equals 1 if X is true, and if X is false. What we are actu- 
ally interested in is the number of solutions that are normalisable for all r 
(Nl(1)). Each such solution will be equal to some combination of the nf 
well behaved solutions to (|8|,||) at large r, and a combination of the N®(1) 
suitable solutions to (10,0]) for small r. If there are only nf, or less, suitable 
small r solutions, then in general any combination of the large r solutions 
will not be well behaved at r = 0. If there are nf + k suitable small r solu- 
tions, then k independent combinations of the large r solutions will be well 
behaved everywhere. It may be possible to get more solutions by fine tuning 
the theory, in which case the index derived would be a lower bound. 
The number of normalisable solutions for a given Z is 



(15) 



where [x] + is defined to be equal to zero if x < 0, and x if x > 0. 

This is not true if the equations obtained from @ can be split into 
several independent sets. This will occur when M a p is a direct product 
of mass matrices. In this case the mass matrix can be split up into smaller 
matrices, which can be analysed individually. Even when M a a is not a direct 
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product of other matrices, it may still be possible to split the equations into 
two independent sets. This case will be considered separately later. 

Since U% and V£ are determined by real equations, each complex so- 
lution gives two real solutions. This suggests that the total number of left 
moving zero modes, Nl, is 2J2i However, as can be seen from (0), 

solutions for I = k and I = —k are equal. For I = U% = ±V^, so one of 
V'a's solutions is zero. Thus the total number of independent real solutions 
is 



E 

i 



£('[■ 



/ < 



.a=l 



+ 1 



I > -Qa + 



- Uf 

+ (16) 

The summation is over all values of I that give single valued ip. Since all the 
Higgs fields which make up M Q/ g are single valued, (||) implies all q Q or all 
q a — 1/2 are integers (assuming M a a is not a product of smaller matrices), 
in which case respectively I — 1/2 or I is an integer. 

A similar analysis can be applied to right moving zero modes. For large 
r the behaviour is the same. For small r, solutions are well behaved if 
I > Qa + 1/2 or / < — q a — 1/2. This gives 



E 

i 



E 



I < 



+ 1 



l > q a + h 



rtf 



(17) 

If q a is positive then one or two of the I [. . .] terms will be non-zero. If q a 
is negative, one or zero of them will be non-zero. By splitting q a into positive 
and negative eigenvalues, and ordering them, (|i"6|) and (|l7]) can be simplified. 
If there are n+ positive and n_ negative q a s then, after reordering, q a = 
(pi, P2 ••• Pn+, -9i, ~92 ■■■ ~9n-, ... 0), where pj,gj > 0. Clearly if 
the string gauge eigenvalues are not integer, there are no zeros. The I [. . .] 
terms in (|l6| ) and (|l7]) can be combined to give 



^ = E 



3=1 



E 1 [-93 ~ 2 < 1 < 9j + \ 

3=1 



n+ 



E 1 + 2 < ; < 9j - h] - E 1 [-Pi - 2 < 1 < PJ + 2 

3=1 3=1 



(18) 



(19) 
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These expressions can be further simplified if the eigenvalues are ordered. 
If pi > p% > ■ ■ ■ > p n+ > and 51 > 52 > • • • > 5n_ > 0, it is possible to 
evaluate the I summation by considering cancellation of the pj and gj terms. 
This gives 

min(n_,n+) n_|_ 

N L = J2 2 \Pj -9j} + + E 2 Pj > ( 20 ) 

j=l j=n-+l 
min(n_ ,n+) n~ 

Nr = E 2 l9j-Pj} + + E 2 ft- • ( 21 ) 
Taking the difference of these results gives 

min(n_,n_|-) n_|_ n_ 

1 A = N L -N R = £ 2 (Pi-9j)+ E 2 ^i- E 2^(22) 

j=l j=n_+l j=n++l 

= E 2 g a = ^[lndetM]to • (23) 

a=l 

This is in agreement with another index theorem obtained elsewhere || [nj . 
The other index theorem was obtained by a different method and only gave 
I A, not Nl and Nr. 

The total number of zero modes is 

min(n_,n+) n+ n_ 

J = jV L + jV R = E 2 \Pj~9j\+ E 2 P; + E 2 9j > (24) 

j=l j=n_+l j=n++l 

where only one of the last 2 terms contributes, depending on whether n + or 
n_ is bigger. This is also true of (22) and ( pO| , p| ) . 

If 2 is to be zero, then for every positive q a there must be one negative 
qp with the same magnitude. If every fermion field couples to a Higgs field 
with winding number zero, this will be the case. 

The above approach fails if C a p is of the form 



J a /3 

K a p 



(25) 



As C a /3 is assumed to have no zero eigenvalues, J a/ 3 and K a p are both n x n 
matrices, where n = nf/2. In this case when (^0) are substituted into ([|), 
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two independent sets of equations are obtained. Expressions for Nl and Nr 
can found by considering just one set of these solutions. Putting 



U^(r)e lW 
V£*(r)e- iW ) 



a 
a 



1 . . . n 

n + 1 . . . nf 



(26) 



and ^) into @ gives (|8],f|) for large r and (10,11) for small r, but with 
a more restricted range on the indices. For large r, n of the rif complex 
solutions are normalisable. For a given I at small r, there is one normalisable 
solution for each q a satisfying I > —q a + 1/2 (a = 1 . . . n) or I < q a — 1/2 
(a = n + 1 . . . rif). 

Matching the solutions for large and small r gives 



i < 



.a=l 



+ 1 



I > -Qa + 



n 



(27) 



The solutions for different / are independent (unlike the previously consid- 
ered cases), so the total number of real solutions (Nl) is just twice the 
number of complex solutions, thus iV^ = 2^ Nl(1)- A similar expression 
can be obtained for Nr. 

The I [. . .] terms in ( |27| ) can be combined. Defining gj = —qj and pj = 
qj+n for j = 1 ... n, (p7|) and the corresponding expression for Nr, become 



Nt 



Nr 



2E 



2E 

l 



n 

Y,{ 1 [9j + l< l < Pj - \] - 1 [Pj - \ < 1 < 9 3 + 5 



71 

E ( r [pj + I < ; < 9j - |] - 1 [gj - \ < I < Pj + 3 



(2?) 
+ 

(29) 



If the string gauge eigenvalues are reordered, so that p\ > ... > p n and 
5i > • • • > 9n, the expressions reduce to (2C) and (21), with n + = n_ = n. 
This is not identical to the previous result, since pj and gj are defined 
differently. 

When there are just two fermion fields involved, all the results reduce to 



I 



%1 + Q2) 



(30) 
(31) 
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3 An 50(10) GUT with Strings 



One example of a phenomologically credible grand unified theory (GUT) 
has the symmetry breaking 



50(10) ^ SU{5) x Z 2 (32) 



$45 
*10 



517(3) C x SU(2) L x U(1) Y x Z 2 (33) 
SU(3) C x U(1) Q x Z 2 . (34) 



The discrete Z 2 symmetry allows the formation of topologically stable cosmic 
strings. One possibility is an abelian string, in which case 

$126 = e in %/(r) , (35) 

where <pQ is the usual VEV of &\ 2 q, /(0) = and /(oo) = 1. To give a 
vanishing covariant derivative at infinity there is also a non-zero gauge field. 
The only stable abelian strings have |n| = 1, but higher winding number 
strings may have long lifetimes. Nonabelian strings also form in this model, 
but they do not have zero modes at high temperatures, and so will be ignored 

i- 

The string gauge field has a non-trivial effect on the electroweak Higgs 
field, $i . Thus 
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H e ime + H' e- ime 



h(r) , (36) 



where Hq and Hq are the usual VEVs of the components of <l?io, h(oo) = 1, 
and h(0) is small (m = 0) or zero (m ^ 0). There is also a non- vanishing 
electroweak gauge field, m is determined by the GUT string, and is equal 
to the nearest integer to — n/5, so the electroweak Higgs field does not wind 
around a topologically stable abelian string. 

Of the three Higgs fields, only $126 an d $10 couple to fermions. Only 
right-handed neutrinos couple to <fo , while neutrinos of either helicity couple 
to Hq. The full mass matrix is a direct product of several matrices. The 
neutrinos get their masses from 



M<» = I m GUT.f(r)e mf? m u h(r)e 
" { m u h(r)e- imd 



-imd 



(37) 



where m u is the up quark mass and ttt-gut ~ |</>o| ~ 10 16 GeV. In the notation 



of (ph ^2 = V L and tpi = ia 2 v^). The gauge eigenvalues of the two fields are 
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then qi = n/2 and q2 = —n/2 — m. The other particles just couple to Hq or 



At high temperatures only $126 is non-zero and m u = 0, so only ipi 
features in the zero mode analysis. Using (|3l| ) gives the number of neutrino 
zero modes as 

Nl = [n] + , N r = [-n} + . (38) 
At low temperatures m u 7^ 0. This time 

N L = [-2m} + , N R = [2m} + . (39) 

Since \2m\ < |n|, some of the zero modes will be destroyed. For a stable 
n = 1 string all zero modes are destroyed. This is in agreement with results 
obtained in ||. Thus, since higher n strings almost certainly decay, there 
are zero modes before, but not after the electroweak phase transition. The 



neutral current in the string disperses [13] and any vortons formed would 
dissipate after about 10 -10 sec [11|. Before the electroweak phase transition 



from about 10 10 GeV - 10 2 GeV the universe would undergo a period of mat- 
ter domination. Once the vortons dissipate there would be some reheating 
of the universe. However the electroweak interactions and physics below the 
phase transition would be unaffected. 

In a more arbitrary U(l)x (Standard Model) theory, the ratio of m and 
n could be greater than 1, in which case extra zero modes would be created 
at the electroweak phase transition. 

The destruction of zero modes at phase transitions will occur in a wide 
range of models. In a theory with majorana mass terms, such as the right 
handed neutrino term above, the first version of the index theorem can be 
used (pi}) . After a phase transition, if every fermion couples to a non- winding 
Higgs field, such as the electroweak Higgs field above, then all previously 
formed zero modes will be destroyed (see remarks after (|24])). 



4 A Model with Persistent Zero Modes 

In this section we consider a model in which zero modes survive a subsequent 
symmetry breaking, despite coupling to a Higgs field with a non-winding 
component. Consider fermions which have the usual Yukawa couplings with 
two Higgs fields. Neglecting gauge fields, the form of the Lagrangian for 
each pair of fermions is then 

Ct = -ifL-fdpfL - ifafdrfn + m(f L 4>*f R + f R <f>f L ) , (40) 
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where <f> is a combination of the two Higgs fields, only one of which gives 
rise to a string. This leads to the Dirac equations 

- h^f L + m<t>*f R = , (41) 

- h^fn + m^f L = . (42) 
Using the following forms for the left and right-handed spinors, 





( 


-A 




<°A 

c 


h = 




f R = 




{ 


-b) 







the equations of motion become 



/ w + k ie- ie (-d r + l -d e 

ie ie (d r + l -d e ) -w + k 

—mcj) 

V —m<b 



We can now try the usual zero mode analysis 
A = D = 0. The equations of motion reduce to 



m<jj* 





\ 







m<j)* 


( A A 


— w + k 


ie- id (d r -^d 


e) 


c 


ie ie {-d T -\d B 


) w + k 


) 





let w 



le 



and 



mo 



d r d e B + m6*C = 







ie ie [d r + -de)C 
r 



(44) 
k and set 



(45) 
(46) 



The Higgs field has two parts, a winding part from the string and a 
constant part from a second symmetry breaking, 



4> = f{r)e ie +p . 



(47) 



Note that this form for <j> cannot be written in the form used in (4) and 
so the index theorem developed in section 2 does not apply in this case. 
Changing variables to X and Y with 



rr 

X= f(p)dp + prcosd + c 
Jo 



(48) 
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(dX i dX 



f{r) +p 



COS I 



db -(— rsinl 

r 



we have 

V dr r 86 

Taking B and C independent of Y the equations of motion become 

- ie~ ie [fir) + p(cos 9 + ism 9))d x B + m(j)*C = 

and 

- m^B - ie i9 [f(r) + p(cos 6 - i sin 6)}d x C = . 
The Higgs fields cancel and we have 

i-iB) tX + mC = , m(-iB) + C tX = 

which have a solution, 



-iS) 
C 



-rnX 



(49) 

(50) 
(51) 

(52) 
(53) 



We have explicitly constructed a zero mode after the second phase transi- 
tion. The X coordinate is similar to the usual radial coordinate, but 'centres' 
on the effective zero of the resultant Higgs field, rather than the core of the 
string. 



5 Index Theorems and Spectral Flow 

We have shown that zero modes can acquire masses at subsequent phase 
transitions. No matter how small this mass the spectrum of the Dirac op- 
erator changes significantly. If we compare the Dirac spectrum with a zero 
mode and a low lying bound state with infinitesimal mass (fig. [I]), we see 
that an arbitrarily small perturbation to the zero mode introduces an entire 
new branch to the spectrum. Any massive state gives a spectrum that is 
symmetric about both the w and k axes, there is always a reference frame 
in which the particle is at rest and others where it is moving up or down the 
string. Conversely the zero mode, which is massless, can only move in one 
direction along the string and its spectrum is asymmetric. The transition 
from zero mode to low lying bound state causes drastic changes in the spec- 
trum and can be brought about by infinitesimal changes in the value of one 
Higgs field. If we consider the species with the zero mode alone, this infinite 
susceptibility to the background fields appears unphysical. However, when 
we include the massless neutrino in the 50(10) model the spectral changes 
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Figure 1: The Dirac spectrum with a zero mode (left) and a very low lying 
bound state (right). Both spectra also have a bound state and continuum. 



are less worrying. For a small coupling between the two neutrinos, both the 
before and after spectra have a continuum of massless or nearly massless 
states. These states can be used to build the extra branch of the perturbed 
zero mode spectrum, allowing small changes in the overall spectrum for 
small changes in the background fields. 

This observation leads us to conjecture that zero modes can be removed 
only if they become mixed with other states. 

The coupling between the left and right handed neutrinos and the elec- 
troweak Higgs field need not be artificially small, the small mass of the light 
neutrino can be generated by the seesaw mechanism [14]. This coupling is 
present prior to the electroweak phase transition and allows transitions of 
the form vj_, + Ur — > //, where / is any light fermion from the standard 
model and the intermediate state is an electroweak Higgs. Such interactions 
allow zero modes on the string to scatter from massless neutrinos in the 
surrounding plasma and provide a current damping mechanism that affects 
current build up prior to the electroweak transition. 
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6 Conclusions 



In this paper we have seen that the microphysics of cosmic strings can be 
influenced by subsequent phase transitions. Fermion zero modes, and con- 
sequently superconductivity, of the strings can be created or destroyed by 
such phase transitions. In determining whether or not a cosmic string is 
superconducting it is not enough to just consider this at formation, but to 
follow the microphysics through the multiple phase transitions that the sys- 
tem undergoes. It is possible for vortons formed at high energy to dissipate 
after a subsequent phase transition if the relevant fermion zero mode does 
not survive the phase transition, thus vortons bounds could be evaded. This 
is the case with vr zero modes on the SO{ 10) string ||. Prior to dissipation 
there could be a period of vorton domination, after the phase transition the 
universe would reheat and then evolve as normal. 

To enable a systematic analysis of this effect we have derived a gener- 
alised index theorem. Our index theorem is especially applicable to theories 
where the fermions acquire mass from more than one Higgs field. We applied 
the index theorem and also considered spectral flow. As a result we conjec- 
ture that zero modes are destroyed when they mix with other fermions that 
acquire mass at a subsequent phase transition from a non-winding Higgs 
field. 
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